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Abstract 

With non-equilibrium applications in mind we present in this paper (first in a 
series of three) a self-contained calculation of the hydrostatic pressure of the 
O(N) A</> 4 theory at finite temperature. By combining the Keldysh-Schwinger 
closed-time path formalism with thermal Dyson-Schwinger equations we com- 
pute in the large N limit the hydrostatic pressure in a fully resumed form. 
We also calculate the high-temperature expansion for the pressure (in D = 4) 
using the Mellin transform technique. The result obtained extends the results 
found by Drummond et al. [1] and Amelino-Camelia and Pi [2]. The latter are 
reproduced in the limits m r (0)^0, T^oo and T^oo, respectively. Impor- 
tant issues of renormalizibility of composite operators at finite temperature 
are addressed and the improved energy-momentum tensor is constructed. The 
utility of the hydrostatic pressure in the non-equilibrium quantum systems is 
discussed. 
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1 Introduction 

In order to give a theoretical description of the properties of matter under extreme 
conditions (like neutron stars, the early universe or heavy-ion collisions) one is of- 
ten forced to use the statistical quantum field theory (QFT). The latter is due to 
inherent quantum nature of these processes and due to an overwhelming number 
of degrees of freedom involved. In recent years, considerable effort has been de- 
voted to the understanding of both equilibrium and non-equilibrium behavior of 
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such systems (see e.g. [3, 4] and citations therein). In fact, the equilibrium descrip- 
tion is worked out relatively well and number of methodologies for doing quantum 
field theory on systems at or near (local) equilibrium is available. On this level 
two modes of description have been formulated: imaginary-time (or Matsubara) 
approach [5, 6, 7, 8] and real-time approach [5, 6, 7, 9]. In contrast to equilib- 
rium, the theoretical understanding of non-equilibrium quantum field theories is 
still very rudimentary. Complications involved are essentially twofold. The first is 
related to an appropriate choice of the non-equilibrium initial-time conditions and 
their implementation into quantum description [4, 10]. The second problem is to 
construct the density matrix pertinent to the level of description one aims at. The 
latter requires usually some sort of coarse-graining (e.g. truncation of higher point 
Wigner functions in the infinite tower of Schwinger-Dyson equations [11]) or pro- 
jecting over irrelevant subsystems (incorporated e.g. via projection operator method 
[12] or maximal entropy - MaXent - prescription [13]). However, when the density 
matrix is known one may, in principle, apply the cummulant expansion to convert 
the calculations into those mimicking usual equilibrium techniques [11, 14]. Yet, the 
boundary problem prohibits per se many of equilibrium approaches. Imaginary-time 
approach is clearly not applicable due to its lack of the explicit time dependence 
and build-in equilibrium (Kubo-Martin-Schwinger) boundary conditions. Among 
the real-time formalisms only the Schwinger-Keldysh or closed-time path formal- 
ism (CTP) [5, 6, 11, 15, 16] and thermo field dynamics (TFD) [5, 17, 18] has found 
a wider utility in non-equilibrium computations. The CTP formulation was con- 
veniently applied, for instance, in the study of non-equilibrium gluon matter [19], 
cosmological back reaction problem [20] or in the time evolution of a non-equilibrium 
chiral phase transition [21]. On the other hand the non-equilibrium TFD was re- 
cently used in deriving the transport equations for dense quantum systems [22], or 
in a study of transport properties of quantum fields with continuous mass spectrum 
[18]. 

To extract information on the underlying field dynamics or on non-equilibrium 
transport characteristics one needs to specify an appropriate set of observables (be 
it conductivity, damping rates, edge temperature jumps, viscosity, etc.). Pressure is 
often one of the key parameters used in the diagnostics of off-equilibrium quantum 
media. Hydrostatic pressure measurements in superfluid He 4 (i.e. in He II phase) 
[23] and in superconductors [24] provide examples. It is thus clear that an extension 
of the pressure calculations to non-equilibrium systems could enhance our predica- 
tive ability in such circcis cis (realistic) phase transitions, early universe cosmology or 
hot fusion dynamics. However, the usual procedure known from equilibrium QFT, 
i.e, calculations based on the partition function or effective potential [7, 25, 26, 27] 
cannot be employed here. This is because the (grand)-canonical potential from 
which the thermodynamic pressure is derived does not exist away from equilibrium. 
Fortunately, more general definition of pressure, not hinging on existence of (grand)- 
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canonical potential, exists. This is the so called hydrostatic pressure and its form 
is deduced from the expectation value of the energy-momentum tensor. It might 
be shown that in thermal equilibrium the (classical) thermodynamic and (classical) 
hydrostatic pressures are identical on account of the (classical) virial theorem [28]. 

In this and two companion papers we aim at clarifying the calculation of the 
hydrostatic pressure away from equilibrium and at studying its bearings to various 
non-equilibrium situations. Calculation of the expectation value of the energy- 
momentum tensor is, however, quite delicate task even in thermal equilibrium as 
computations involved are qualitatively very different from those known, for in- 
stance, from the effective action approach. This is because the energy-momentum 
tensor is a composite operator and as such it requires a different methodology of 
treatment including a different approach to renormalization issues [5, 29]. It should 
then come as a no surprise that in thermal QFT the equivalence between hydrostatic 
and thermodynamic pressure (or effective action) is more fragile than in correspond- 
ing classical statistical systems. In fact, the validity of the quantum virial theorem 
is by no means established conclusively, and it is conjectured that it could break 
down, for instance, in gauge theories [5]. Besides, there is clearly no virial theorem 
away from equilibrium (not even classically) and so in such a case one must expect 
disparity between hydrostatic pressure and effective action. 

In order to understand the difficulties involved we concentrate in the present 
paper on the calculation of the hydrostatic pressure in thermal equilibrium. To 
this end, we utilize the CTP approach which both in spirit and in many technical 
details mimics the realistic non-equilibrium calculations [11, 13, 21, 14]. Presented 
CTP formalism in addition to its theoretical structure which is interesting in its 
own right, is important because it can be with a minor changes directly applied to 
translationally invariant non-equilibrium QFT systems [13]. In order to keep the 
discussion as simple as possible we illustrate our reasonings on O(N) symmetric 
scalar A 4 theory. The model is sufficiently simple yet complex enough to serve as 
an illustration of basic characteristics of the presented method in contrast to other 
ones in use. The latter has the undeniable merit of being exactly solvable in the 
large- N limit both at zero and finite temperature [25, 30, 31, 32, 33, 34, 35]. It 
might be shown that the leading order approximation in 1/N is closely related to 
the Hartree-Fock mean field approximation which has been much studied in nuclear, 
many-body, atomic and molecular chemistry applications [25, 36]. In addition, in 
the case of a pure state it corresponds to a Gaussian ansatz for the Schrodinger wave 
functional [37]. We will amplify some of these points in later papers. We should also 
emphasize that although the O(N) 4 theory frequently serves as a useful playground 
for study of finite-temperature phase transitions with a scalar order parameter, this 
point is not objective of this work and hence we will not pursue it here. 

The set-up of the paper is the following: In Sections 2 we briefly review the 
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derivation of the thermodynamic and hydrostatic pressures. In Section 3 we lay 
down the mathematical framework needed for the finite-temperature renormaliza- 
tion of the energy-momentum tensor (for an extensive review on renormalization 
of composite operators the reader may consult e.g., refs. [5, 38, 39]). The latter is 
discussed on the O(N) 4 theory. It is a common wisdom that the zero temper- 
ature renormalization takes care also of the UV divergences of the corresponding 
finite temperature theory [5, 7, 40]. The situation with energy-momentum tensor 
is, however, more complicated as there is no well defined expectation value of the 
stress tensor at T = [5, 29]. We show how this problem can be amended at fi- 
nite temperature. The key original results obtained here is the prescription for the 
improved energy-momentum tensor of the O(N) 4 theory. The latter is achieved 
by means of the Zimmerman forest formula. With the help of the improved stress 
tensor we are able to find the corresponding QFT extension of hydrostatic pressure 
and hence obtain the prescription for the renormalized pressure. This latter result 
is also original finding. As a byproduct we renormalize (ft 2 and (ft a <ftb operators. 

Resumed form for the pressure in the large- N limit, together with the discus- 
sion of both coupling constant and mass renormalization is presented in Section 
4. Calculations are substantially simplified by use of the thermal Dyson-Schwinger 
equations. For simplicity's sake our analysis is confined to the part of the parameter 
space where the ground state at large N has the O(N) symmetry of the original 
Lagrangian and the spontaneous symmetry breakdown and Goldstone phenomena 
are not possible (Bardeen and Moshe's parameter space [33]). 

In Section 5 we end up with the high-temperature expansion of the pressure. 
Calculations are performed in D = 4 both for massive and massless fields, and the 
result is expressed in terms of the renormalized mass m r (T) and the thermal mass 
shift 5m 2 (T). The expansion is done by means of the Mellin transform technique. 
In appropriate limits we recover the results of Drummond et al. [1] and Amelino- 
Camelia and Pi [2] for thermodynamic pressure (effective action). 

The paper is furnished with two appendices. In Appendix A we clarify some 
mathematical manipulations needed in Section 4. For the completeness' sake we 
compute in Appendix B the high-temperature expansion of the thermal-mass shift 
Sm 2 (T) which will prove useful in Section 5. 

2 Hydrostatic pressure 

In thermal quantum field theory where one deals with systems in thermal equilibrium 
there is an easy prescription for a pressure calculation. The latter is based on the 
observation that for thermally equilibrated systems the grand canonical partition 
function Z is given as 

Z = e-P n = Tr{e~ pi - H ~^ Ni) ) , (2.1) 
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where f2 is the grand canonical potential, H is the Hamiltonian, Ni are conserved 
charges, are corresponding chemical potentials, and (3 is the inverse temperature: 
P = 1/T (k B = 1). Using identity f3-§^ = -T£> together with (2.1) one gets 

= n - E + fuNi , (2.2) 




with E and V being the averaged energy and volume of the system respectively. A 
comparison of (2.2) with a corresponding thermodynamic expression for the grand 
canonical potential [5, 41, 42, 43] requires that entropy S — — ^, so that 



dVl = -SdT - pdV - Nidfa p=-l— \ . (2.3) 

For large systems one can usually neglect surface effects so E and N { become exten- 
sive quantities. Eq.(2.1) then immediately implies that fl is extensive quantity as 
well, so (2.3) simplifies to 

p --v = W (2A) 

The pressure defined by Eq.(2.4) is so called thermodynamic pressure. 

Since InZ can be systematically calculated summing up all connected closed 
diagrams (i.e. bubble diagrams) [5, 44, 45], the pressure calculated via (2.4) enjoys 
a considerable popularity [1, 5, 6, 46]. Unfortunately, the latter procedure can not 
be extended to out of equilibrium as there is, in general, no definition of the partition 
function Z nor grand-canonical potential Q away from an equilibrium. 

Yet another, alternative definition of a pressure not hinging on thermodynamics 
can be provided; namely the hydrostatic pressure which is formulated through the 
energy-momentum tensor Miy . The formal argument leading to the hydrostatic 
pressure in D space-time dimensions is based on the observation that (O^rr)) is 
the mean (or macroscopic) density of momenta p- 7 in the point x^ 1 . Let P be the 
mean total (D — l)-momentum of an infinitesimal volume V^ ^ 1 ^ centered at x, 
then the rate of change of j-component of P reads 



dt 



[ d D -^ A(e V,x')>= E / ds l (Q^). (2.5) 



In the second equality we have exploited the continuity equation for (0 W ') and 
successively we have used Gauss's theorem 1 . The dV^ D ^ corresponds to the surface 
ofl/^ 1 ). 



1 The macroscopic conservation law for (9'"') (i.e. the continuity equation) has to be postulated. 
For some systems, however, the later can be directly derived from the corresponding microscopic 
conservation law [47]. 
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Anticipating a system out of equilibrium, we must assume a non-trivial distri- 
bution of the mean particle four-velocity U^{x) (hydrodynamic velocity). Now, a 
pressure is by definition a scalar quantity. This particularly means that it should not 
depend on the hydrodynamic velocity. We must thus go to the local rest frame and 
evaluate pressure there. However, in the local rest frame, unlike the equilibrium, 
the notion of a pressure acting equally in all directions is lost. In order to retain the 
scalar character of pressure, one customarily defines the pressure at a point (in the 
following denoted as p(x)) [48], which is simply the "averaged pressure" 2 over all 
directions at a given point. In the local rest frame Eq.(2.5) describes j-component 
of the force exerted by the medium on the infinitesimal volume V^ D ~ l \ (By defi- 
nition there is no contribution to caused by the particle convection through 
.) Averaging the LHS of (2.5) over all directions of the normal n(x), we get 3 

1 D ~ l r 

= Y / ds (G l ,(x')) , (2.6) 

where dQ(n) is an element of solid angle about n and Si~ 2 is the surface of (D — 2)- 
sphere with unit radius (/ dQ(n) = Si~ 2 = 2ir~5~ /r(-^l)) . On the other hand, 
from the definition of the pressure at a point x^ we might write 

(srr % s ^ * - -«*) u„ * • w 

here the minus sign reflects that the force responsible for a compression (conven- 
tionally assigned as a positive pressure) has reversed orientation than the surface 
normals n (pointing outward). In order to keep track with the standard text-book 
definition of a sign of a pressure [42, 48] we have used in (2.7) the normal n in a 
contravariant notation (note, n* = — iij). Comparing (2.6) with (2.7) we can write 
for a sufficiently small volume V^ ^ 

p(*) = -p^£W0)- (2-8) 

We should point out that in equilibrium the thermodynamic pressure is usually 
identified with the hydrostatic one via the virial theorem [5, 49]. In the remainder 
of this note we shall deal with the hydrostatic pressure at equilibrium. We shall 
denote the foregoing as V(T), where T stands for temperature. We consider the 
non-equilibrium case in a future paper. 

2 To be precise, we should talk about averaging the normal components of stress [48]. 

3 The angular average is standardly defined for scalars (say, A) as; / A df2(n)/ J dCl(n), and for 
vectors (say, A 1 ) as; J2j I ^ nJ <^( n )/ / dQ(n). Similarly we might write the angular averages 
for tensors of a higher rank. 
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3 Renormalization 



If we proceed with (2.8) to QFT this leads to the notorious difficulties connected 
with the fact that Miy is a (local) composite operator. If only a free theory would be 
in question then the normal ordering prescription would be sufficient to render (G M ^) 
finite. In the general case, when the interacting theory is of interest, one must work 
with the Zimmerman normal ordering prescription instead. Let us demonstrate the 
latter on the O(N) <fi 4 theory. (In this Section we keep N arbitrary.) Such a theory 
is defined by the bare Lagrange function 

£> = \ £ ((^) 2 - "M) - ^ [Ewj , (3-i) 

and we assume that wig > 0. The corresponding canonical energy-momentum tensor 
is given by 

©r = £^ a d> a -<r>c. (3.2) 

a 

The Feynman rules for Green's functions with the energy-momentum insertion can 
be easily explained in momentum space. In the reasonings to follow we shall need 
the (thermal) composite Green's function 4 

D^(x n \y) = (T* {^(xx) . . • 0r(x n )er (y)}> . (3.3) 

Here the subscript r denotes the renormalized fields in the Heisenberg picture (the 
internal indices are suppressed) and T* is the so called T* product (or covariant T 
product) [26, 50, 51, 52]. The T* product is defined in such a way that it is simply the 
T product with all differential operators V^. pulled out of the T-ordering symbol, 
i.e. 

T*{V*\M*l) ■ -.V^MXn)} = V(id M )T{<j> r (xi) ■ ..<t>r{Xn)} , (3.4) 

where V{id^) is just a useful short-hand notation for V X ^V X £ . . . ZF™. In the case 

of thermal Green's functions, T represents a contour ordering symbol [5, 6, 7]. It is 
the mean value of the T* ordered fields rather than the T ones, which corresponds 
at T = and at equilibrium to the Feynman path integral representation of Green's 
functions [52, 53]. 

A typical contribution to Q% v (y) can be written as 

P m 0(y)^ 2 0(y)...P M? >(y), (3.5) 
so the typical term in (3.3) is 

V(td M ) (T'{<f> r (xi)---Mx n)<Kvi) ■ -.<KVk)}) \ yi=y = D M (x n \y k ) \ m=y . 

4 By <j> we shall mean the field in the Heisenberg picture. The subscript H will be introduced in 
cases when a possible ambiguity could occur. 
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Performing the Fourier transform in (3.3) we get 



D^(p n \p)= E [(n^)fr) D 6 D (p-j: qj )V% ) ( q{lt} ) D(p"\q k ) , (3.6) 

fc={2,4} J \i=l y Z7l > ) j=l 



where T>^(. . .) is a Fourier transformed differential operator corresponding to the 
quadratic (k = 2) and quartic (k = 4) terms in 0£". Denoting the new vertex 
corresponding to T>^(. . .) as ®, we can graphically represent (3.3) through (3.6) as 



D^(p n \p) = 





Figure 1: The graphical representation of D tiy (p n \p). 



For the case at hand one can easily read off from (3.2) an explicit form of the bare 
composite vertices, the foregoing are 



jxv ~ V^(q M ) = \ 5 ab {2( qi -p)»cK - g^((qi -p)xql - m 2 )} 
q i a / C> \q 2 



P 



d ~ ^(4) = 1 sW L {gabbed + $aJbd + ^ad^bc) ~ 5S ab 5 cd 5 ac } 



(For the internal indices we do not adopt Einstein's summation convention.) The 
blobs in Fig. 1 comprise the sum of all n + 2- and n + 4- (not necessarily connected) 
Green functions. As usual, the disjoint bubble diagrams in Green functions (blobs) 
can been divided out from the very beginning. We have also implicitly assumed that 
the summation over internal indices is understood. 
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In case when we deal with finite temperature, we choose the contour ordering in 
(3.3) to run along the time contour depicted in Fig.3. It is possible to show that 
for Green's function calculations only horizontal paths contribute [16, 54, 55]. In 
addition, the "physical" fields occurring on the external lines of Green's functions 
have time arguments on the upper horizontal path (type-1 fields) while the "ghost" 
fields have time arguments on the lower horizontal path (type-2 fields). The latter 
modify the Feynman rules in a nontrivial fashion [6, 7, 16]. From the foregoing 
discussion should be clear that in the case of thermal composite Green's function, 
the new (composite) vertices are of type-1 as the fields from which they are deduced 
are all physical 5 . 

3.1 Renormalization of (f) a (x)(f)b(x) 

Now, if there would be no insertion in (3.3), the latter would be finite, and so it is 
natural to define the renormalized energy-momentum tensor [0^1 ( or Zimmermann 
normal ordering) in such a way that 



is finite for any n > 0. To see what is involved, we illustrate the mechanism of 
the composite operator renormalization on <f) a (x)4>b{x) first - the energy-momentum 
tensor case will be postponed to Section 3.2. In the following we shall use the mass- 
independent renormalization, and for definiteness we chose the minimal subtraction 
scheme (MS). In MS we can expand the bare parameters into the Laurent series 
which has a simple form [26, 39, 53], namely 



Here ao and bo are analytic in D = 4. The parameter [i is the scale introduced by 
the renormalization in order to keep A r dimensionless. An important point is that 
both cifc's and fr^'s are mass, temperature and momentum independent. 

It was Zimmermann who first realized that the forest formula known from the 
ordinary Green's function renormalization [26, 38] can be also utilized for the com- 
posite Green's functions rendering them finite [38, 56]. That is, we start with Feyn- 
man diagrams expressed in terms of physical (i.e. finite) coupling constants and 
masses. As we calculate diagrams to a given order, we meet UV divergences which 
might be cancelled by adding counterterm diagrams. The forest formula then pre- 
scribes how to systematically cancel all the UV loop divergences by counterterms 

5 For a brief introduction to the real-time formalism in thermal QFT see for example [5, 6, 9]. 



Dr(x n \y) = (T* {<j> r (xi) ■ ■ ■ <PrM [6f ]}> , 




(3.7) 




(3.8) 
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to all orders. However, in contrast to the coupling constant renormalization, the 
composite vertex need not to be renormalized multiplicatively. We shall illustrate 
this fact in the sequel. Let us also observe that in the lowest order (no loop) the 
renormalized composite vertex equals to the bare one, and so to that order A = [A], 
for any composite operator A. 

Now, from (3.7) and (3.8) follows that for any function F = F(m r , A r ) we have 



dF dml dF m 2 dF 



dmj dw? r dm\ m 2 dm^ 

So particularly for 

F — D(xi, . . . , X n ) — (T*{0 r (xi) . . . (frr(Xn)}) , 

one reads 



7z- 



Jd D xJ2 jv<t>^r{xi)...<t>r{Xn) ^(x) exp(^[0, T]) 
N 

J d D x D a(x 1 ,...,x n \x;ml). (3.9) 



Here A/" -1 is the standard denominator of the path integral representation of Green's 
function. We should apply the derivative also on M but this would produce discon- 
nected graphs with bubble diagrams. The former precisely cancel the very same 
disconnected graphs in the first term, so we are finally left with no bubble diagrams 
in (3.9). In the Fourier space (3.9) reads 



/ i\ N 

D( Pl ,...,p n )= (-2JE D a ( Pl ,...,p n \0;ml). (3.10) 



m 2_d_ ' : A 

r dml 



As the LHS is finite there cannot be any pole terms on the RHS either, and so 
X)a m o^a i s by itself a renormalized composite operator. We see that m\ precisely 
compensates the singularity of X^Li 0a- 

Now, it is well known that any second-rank tensor (say M ao ) can be generally 
decomposed into three irreducible tensors; an antisymmetric tensor, a symmetric 
traceless tensor and an invariant tensor. Let us set M ab = 4> a 4>b, so the symmetric 
traceless tensor K ab reads 

n 

K ab (x) = Mx)M%) - Sab/N ]>>c(*) , (3.11) 

c=l 
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whilst the invariant tensor I ab is 



N 

I ab {x) = 5 ab /NY J 4>l{x). 

c=l 

Because the renormalized composite operators have to preserve a tensorial structure 
of the bare ones, we immediately have that 

K ab = A 1 [K ab ] and I ab = A 2 [I ab \, (3.12) 

where both A\ and A2 must have structure (1 + XI (poles)). The foregoing guarantees 
that to the lowest order K a b = [K ab ] and l a b = [lab]- As we saw in (3.10), m^I ao is 
renormalized, and so from (3.12) follows that mll ab = C [I a b]- Here C has dimension 
[m 2 ] and is analytic in D = 4. We can uniquely set C = m 2 because only this choice 
fulfils the lowest order condition I ao = [I ao ] (c.f. Eq.(3.8)). Collecting our results 
together we might write 



c 



E 



2 

c 
c 



^E^c], (3.13) 



C 



2 

with Zz<p2 = A 2 = ~^pr ' In the second equality we have used an obvious linearity 
[38] of [...]. From (3.11) and (3.13) follows that 

5 N 

Mx)M*) = AMaWM*)] - -tf( a i - YMW\ ■ ( 3 - 14 ) 

iV c=l 

So particularly for (p 2 one reads 

<Pl = ^ ((N - 1)A, + Z^) [d>l} -jf(A 1 - Z^) E[0 C 2 ] • (3.15) 

From the discussion above it does not seem to be possible to obtain more information 
about A\ without doing an explicit perturbative calculations, however it is easy to 
demonstrate that A\ 7^ Z^. To show this, let us consider the simplest non-trivial 
case; i.e. N=2, and calculate A 1 to order A r . For that we need to discuss the 
renormalization of the n-point composite Green's function with, say, <f> 2 insertion. 
To do that, it suffices to discuss the renormalization of the corresponding 1PI n-point 
Green's function. The perturbative expansion for the composite vertex to order A r 
can be easily generated via the Dyson-Schwinger (DS) equation [57] and it reads 6 

6 Throughout the paper we accept the usual convention: Ordinary (not necessarily connected) 
N-point Green's functions are represented with dotted blobs with N external legs, connected N- 
point Green's functions are represented with hatched blobs with N external legs and 1-PI N-point 
Green's functions are represented with cross hatched blobs with N truncated legs (represented by 
solid circles in vertices). 
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Here cross-hatched blobs refer to (renormalized) 1PI (n + 2)-point Green's func- 
tion, circled indices mark a type of the field propagated on the indicated line, and 
uncircled numbers refer to thermal indices (we explicitly indicate only relevant ther- 
mal indices). The counterterms, symbolized by a heavy dot, are extracted from 
the boxed diagrams (elementary Zimmermann forests). In MS scheme one gets the 
following results: 

4^ / (£p {Dll(g)Dll(-g) " B 12 (g)D 12 (-g)} | MS pole term 

\d ml (^#1 A r m?~ 2 ) | MS = -\rfi 4 - D /2 (D - 4) (4vr) 2 
X r /j 4 - D /Q (D - 4) (6tt) 2 . 



Here E>n and E>i 2 are the usual thermal propagators in the real-time formalism 
[5, 6, 9] (see also Section 4). From (3.16) we can directly read off that 

As the coefficient before 2 * s n °t zer o, we conclude that A\ 7^ It is not a 

great challenge to repeat the previous calculations for the 0i0 2 insertion. The latter 
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gives 



X n 4 " D 



3 (D - 4) (4tt) 2 

Eq.(3.15) exhibits the so called operator mixing [26]; the renormalization of (j) 2 a 
cannot be considered independently of the renormalization of <fi 2 (c ^ a). The latter 
is a general feature of composite operator renormalization. Note, however, that 4> a 4>b 
(a 7^ b) do not mix by renormalization, i.e. they renormalize multiplicatively. It 
can be shown that composite operators mix under renormalization only with those 
composite operators which have dimension less or equal [26, 38, 56]. 

Unfortunately, if we apply the previous arguments to n — 0, the result is not 
finite; another additional renormalization must be performed. The fact that the 
expectation values of [. . .] are generally UV divergent, in spite of being finite for the 
composite Green's functions 7 , can be nicely illustrated with the composite operator 
[(f) 2 ] in the N — 1 theory. Taking the diagrams for -D(0|0) and applying successively 
the (unrenormalized) DS equation [5, 57] we get 




(3.17) 



Eq.(3.17) might be rewritten as 
£>(0|0) = £(010)1,0 

+ TQ I \^ 5D{ % qj) DamP ^\* D(q 6 ), (3.18) 

where D amp (q m \0)\ X k is the m-point amputated composite Green's function to order 
Ajr, and D(q m ) is the full m-point Green's function. The crucial point is that we 
can write D(0|0) as a sum of terms, which, apart from the first (free field) diagram, 
are factorized to the product of the composite Green's function with n > and 
the full Green's function. (The factorization is represented in (3.17) by the dashed 
lines.) Note that the expansion (3.17) is not unique as various other ways of pulling 

7 Also called the matrix elements of [...]. 
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vertices out of Green's function may be utilized but this particular form will prove 
to be important in the next section (see Eq.(3.26)). 

Now, utilizing the counterterm renormalization to the last two diagrams in (3.17) 
we get situation depicted in Fig.2. Terms inside of the parentheses are finite, this 
is because both the composite Green's functions (n > 2 !) and the full Green's 
functions are finite after renormalization. The counterterm diagrams, which appear 
on the RHS of the parentheses, precisely cancel the UV divergences coming from the 
loop integrations over momenta q\ . . . qi which must be finally performed. The heavy 
dots schematically indicates the corresponding counterterms. In the spirit of the 
counterterm renormalization we should finally subtract the counterterm associated 
with the overall superficial divergence 8 related to the diagrams in question. But as 
we saw this is not necessary; individual counterterm diagrams (Zimmermann forests) 
mutually cancel their divergences leaving behind a finite result. 




Figure 2: Counterterm renormalization of the last two diagrams in Eq.(3.17). (Cut 
legs indicate amputations.) 



So the only UV divergence in Eq.(3.17) which cannot be cured by existing coun- 
terterms is that coming from the first (i.e. free field or ring) diagram. The foregoing 
divergence is evidently temperature independent (to see that, simply use an explicit 
form of the free thermal propagator E>n). Hence, if we define 

(0 2 )reno rm = ([<t> 2 ]} ~ (O|[0 2 ]|O), (3.19) 

8 A simple power counting in the (j) A theory reveals [26] that for a composite operator A with 
dimension uia the superficial degree of divergence w corresponding to an n-point diagram is w = 
to a — n. 
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or, alternatively 

(0 2 )renorm = ([0 2 ]) - ([0 2 ])|free fields , ( 3 - 20 ) 
we get finite quantities, as desired. On the other hand, we should emphasize that 

(0 2 > - (O|0 2 |O> = Z e {([0 2 ]> - (O|[0 2 ]|O)} ^ finite in D=A. (3.21) 

An extension of the previous reasonings to any iV > 1 is straightforward, only 
difference is that we must deal with operator mixing which makes (3.19) and (3.20) 
less trivial. 

The important lesson which we have learnt here is that the naive "double dotted" 
normal product (i.e. subtraction of the vacuum expectation value from a given 
operator) does not generally give a finite result. The former is perfectly suited 
for the free theory = 1) but in the interacting case we must resort to the 

prescription (3.19) or (3.20) instead. 



3.2 Renormalization of the energy— momentum tensor 

In order to calculate the hydrostatic pressure, we need to find such (0^)|renorm 
which apart from being finite is also consistent with our derivation of the hydrostatic 
pressure introduced in the introductory Section. In view of the previous treatment, 
we however cannot, however, expect that 0£" will be renormalized multiplicatively. 
Instead, new terms with a different structure than 0{f" itself will be generated during 
renormalization. The latter must add up to 0£" in order to render D^ u (x n \y) finite 9 . 

Now, the key ingredient exploited in Eq.(2.5) is the conservation law (continuity 
equation). It is well known that one can "modify" 0^ in such a way that the new 
tensor Q^ u preserves the convergence properties of 0{f". Such a modification (the 
Pauli transformation) reads 

9"" = 9f + d x X x ^ , 

For scalar fields (3.22) is the only transformation which neither changes the diver- 
gence properties of nor the generators of the Poincare group constructed out 
of 0£" [5, 26, 47, 51]. Because the renormalized (or improved) energy momentum 
tensor must be conserved (otherwise theory would be anomalous), it has to mix with 
0£" under renormalization only via the Pauli transformation, i.e. 

[9f] = 0f + d x X x » u . (3.23) 

9 In fact it can be shown [5, 38] that the Noether currents corresponding to a given internal 
symmetry are renormalized, i.e. J a — [J a ], however, this is not the case for the Noether currents 
corresponding to external symmetries (like 0£" is). 
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In order to determine X x ^, we should realize that its role is to cancel divergences 
present in 0£". Such a cancellation can be, however, performed only by means of 
composite operators which are even in the number of fields (note that is even 
in fields and Green's functions with the odd number of fields vanish). Recalling the 
condition that renormalization can mix only operators with dimension less or equal, 
we see that the dimension of X Xilv must be D — 1, and that X XiiV must be quadratic 
in fields. The only possible form which is compatible with tensorial structure (3.22) 
is then 

N 

X X ^= ]T c ab (\ r ;D) (dV-dV) Mb- (3.24) 

a,b=l 

From the fact that Q£ u and [Q£ u ] are O(N) invariant (see Eq.(3.2)), d x X x ^ u must 
be also O(N) invariant, so c a b = 8 a bC- Thus, finally we can write 

N 

[©H = ©r + 4K; D) ]T (>cT - g^d 2 ) <Pl , (3.25) 

a=l 

with c = Co + XXpoles); here Co is analytic in D. Structure of c(A r ; D) could be 
further determined, similarly as in the N — 1 theory, employing a renormalization 
group equation [39]. We do not intend to do that as the detailed structure of c 
will show totally irrelevant for the following discussion, however, it turns out to be 
important in non-equilibrium cases. 

Now, similarly as before, [0^] gives the finite composite Green's functions if 
n > but the expectation value ([0^]) is divergent (for discussion of the N — 1 
theory see, e.g., Brown [39]). The unrenormalized DS equation for D^(0\0) reads [5] 




The structure of the composite vertices in (3.26) is that described at the beginning 
of this Section. Note that the amputated composite Green's functions in individual 
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parentheses are of the same order in A r . Performing the counterterm renormaliza- 
tion as in the case of ([0 2 ]), we factorize the graphs inside of parentheses into the 
product of the renormalized 2- (and 6-) point composite Green's function and the 
renormalized full 2- (and 6-) point Green's function. The latter are finite. The UV 
divergences arisen during the integrations over momenta connecting both compos- 
ite and full Green's functions are precisely cancelled by the remaining counterterm 
diagrams. Only divergence comes from the free-field contribution, more precisely 
from the T = ring diagram. Defining 

(©D|renor m =([er])-(0|[er]|0), (3.27) 

or 

(©Dlrenorm = ([©f ]> - ([©f ]>lfree field , (3-28) 

we get the finite expressions. Note that the conservation law is manifest in both 
cases. In equilibrium (and in T = 0) we can, due to space-time translational 
invariance of (...), write 

([©f ]> = (©f > + d x (X x ^) = (0f ) . (3.29) 

Using (3.27) or (3.28) we get either the thermal interaction pressure or the interac- 
tion pressure, respectively. This can be explicitly written as 

nh.int.Cn = V(T) -V{0) = -p^T) E{< e ci> - (°l©cJ0)} , (3.30) 

or 

= V{T) - V bcc field (T) = - p3T) E {< "> " (©c^lfree field} • (3-31) 

In order to keep connection with calculations done by Drummond et al. in [1] 
we shall in the sequel deal with the thermal interaction pressure only. If instead 
of an equilibrium, a non-equilibrium medium would be in question, translational 
invariance of (. . .) might be lost, in that case either prescription (3.27) or (3.28) is 
obligatory, and consequently c(A r ; D) in (3.25) must be further specified. 

4 Hydrostatic pressure - calculation 

In the previous section we have prepared ground for a hydrostatic pressure calcula- 
tions. In this section we aim to apply the previous results to the massive O(N) 4 
theory in the large-iV limit. Anticipating an out of equilibrium application, we 
shall use the real-time formalism even if the imaginary-time one is more natural 
in the equilibrium context. As we aim to evaluate the hydrostatic pressure in 4 
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dimensions, we use here, similarly as in the previous Section, the usual dimensional 
regularization to regulate the theory (i.e. here and throughout we keep D slightly 
away from the physical value D — 4). 

In order to actually pursue the pressure calculation we feel it is necessary to 
briefly review the mass and coupling renormalization of the model at hand. This 
will also help to clarify the notation used. While we hope to provide all essentials 
requisite for our task, good discussion of alternative approaches and renormalization 
prescriptions may be obtained for instance in [31, 32, 35]. 

4.1 Mass renormalization 

In the Dyson multiplicative renormalization the fact that the complete propagator 
has a pole at the physical mass leads to the usual mass renormalization prescription 
[26]: " 

m 2 r = m 2 + E(m 2 r ) , (4.1) 

where m r is renormalized mass and S(m^) is the proper self-energy evaluated at 
the mass shell; p 2 = m 2 . In fact, Eq.(4.1) is nothing but the statement that 2-point 
vertex function evaluated at the mass-shell must vanish. The Dyson-Schwinger 
equation corresponding to the proper self-energy reads [31, 40, 57, 58]: where 




S ac |^ c = ; X - xjn 

hatched blobs represent 2-point connected Green's functions whilst cross-hatched 
blobs represent proper vertices I^ 4 ) (i.e. 1PI 4-point Green's function). As S aa are 
the same for all a, we shall simplify notation and write E instead. In the sequel the 
following convention is accepted: 



d 
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The second term in (4.2) actually does not contribute in the large- N limit. It is easy 
to see that the third term does not contribute either. This is because each hatched 
blob behaves at most as iV° whilst T^ 4 ) goes maximally 10 as N^ 1 . Consequently, 
various contributions from the first graph in (4.2) contribute at most N°, whereas in 
the second graph the contributions contribute up to order iV -1 . So the first diagram 
dominates, provided we retain only such 2-point connected Green's functions which 
are proportional to iV° (as mentioned in the footnote, these are comprised only of 
tadpole loops.). After neglecting the "setting sun" graph, Eg. (4.2) generates upon 
iterating the so called superdaisy diagrams [1, 27, 40]. 

Let us now define T,(m 2 ) = A M(m 2 r ). Because the tadpole diagram in (4.2) 
can be easily resumed we observe that 

2\ 1 f d D q i If d D q i 

M{m r ) = -J g2 _ m 2_ S ( m 2) +lt = 2J j2n)D q 2_ m 2 + ie ' ( - 3 ) 

hence we see that S is external-momentum independent. If we had started with 
the renormalization prescription: iTf\p 2 = 0) = — m 2 , we would arrived at (4.1) as 
well (this is not the case for N = 1!). 

At finite temperature the strategy is analogous. Due to a doubling of degrees of 
freedom, the full propagator is a 2 x 2 matrix. The latter satisfies, similarly as at 
T = 0, Dyson's equation 

D = D F + D F (-iS)D. (4.4) 

An important point is that there exists a real, non-singular matrix M (Bogoliubov 
matrix) [5, 6, 9] having a property that 

E - = m ( 1 « F -Ly ) M md E = M " ( t 4 ) M " ■ (45) 

Here Ap is the standard Feynman propagator and * denotes the complex conjuga- 
tion. Consequently, the full matrix propagator may be written as 




D = M p -™o--T+*e _. M . (4.6) 



Similarly as in many body systems, the position of the (real) pole of ID in p 2 fixes 
the temperature-dependent effective mass m r (T) [6, 59]. The latter is determined 
by the equation 

m 2 r {T) = ml + Re (s T (m^(T))) . (4.7) 

10 In the </> 4 theory there is a simple relation between the number of loops (L), vertices (V) 
and external lines (E); 4V = 21 + E. Together with the Euler relation for connected graphs; 
L = I — V + 1 (here I is the number of internal lines), we have L — V = As each loop 

carries maximally a factor of N (this is saturated only for "tadpole" loops) and each vertex carries 
a factor of iV _1 , the overall blob contribution behaves at most as N L ~ V = N~ ~ . 
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From the explicit form of M it is possible to show [5, 6] that ReSn = ReS T . As 
before, the structure of the proper self-energy can be deduced from the correspond- 
ing Dyson-Schwinger equation. Following the usual real-time formalism convention 
(type-1 vertex ~ — i\ , type-2 vertex ~ i\ ), the former reads: 



-iSn = 4 



where 





(4.8) 





ill (yy. \ i i 12 r<xA 2 l 
+ 






1 1 K'/9\l 1 K /X\ 1 1 
+ 4>>// ± #Y/ • + . . . 



and similarly for ID 22- In (4.8) we have omitted diagrams which are of order 0(1/ N) 
or less. Note that the fact that no setting sun diagrams are present implies that the 
off-diagonal elements of S are zero. Inspection of Eq.(4.8) reveals that 

dP q , . . _ A f d D q 



^11 — — 



/p^D llte r) and E 22 = -^/^LD 22 (,;T). (4.9) 

It directly follows from Eq.(4.9) that both S n and S 2 2 are external-momentum 
independent and real 11 . If we define S T (m^(T)) = A M.T(m%(T)), then Eq.(4.7) 
through Eq.(4.9) implies that 

ml(T) = m 2 + A M T (m 2 r (T)) . (4.10) 

A resumed version of Du is easily obtainable from (4.6) [5, 6] and consequently 
(4.9) yields 



d q s(q ) 



Im 



(2%) D C?o/3 - 1 g 2 _ m *(T) + ie 




11 Reality of En can be most easily seen from the largest-time equation [58]. The LTE states 
that £n + S 2 2 + S12 + S 2 i = 0. Because no setting sun graphs are present, X12 + £21 = 0, on 
the other hand En + S 2 2 = 2iImXn (see (4.5)). 
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Let us remark that (4.11) is manifestly independent of any particular real-time 
formalism version. 

In passing it may be mentioned that because Sn(mJ) is momentum independent, 
the wave function renormalization — 1. (The Kallen-Lehmann representation 
requires the renormalized propagator to have a pole of residue i at p 2 = m 2 . The 
former in turn implies that = (1 — S' 11 (p 2 )| p 2 =m 2) _1 = 1.) Trivial consequence of 
the foregoing fact is that if) = and = 

4.2 Coupling constant renormalization 

Let us choose the coupling constant to be defined at T = 0. This will have the 
advantage that the high temperature expansion of the pressure (see Section 5) will 
become more transparent. In addition, such a choice will allow us to select safely 
the part of the parameter space in which spontaneous symmetry breakdown is not 
possible. An alternative renormalization procedure based on the affective action is 
presented in [31]. 

By assumption the fields 4> a have non-vanishing masses, so we can safely choose 
the renormalization prescription for A r at s = (s is the standard Mandelstam 
variable). For example, one may require that for the scattering aa — > bb 

(s = 0) = -X r /N, (b^a). (4. 12) 

The formula (4.12) clearly agrees with the tree level value T^ e ^ abh {s = 0) = —Xq/N. 
Let us also mention that Ward's identities corresponding to the internal O(N) sym- 
metry enforce aaaa to obey the constraint 12 

r^ aaaa ( Pl ;p 2 ;p 3 -p 4 ) = T^ bbaa (p 1 ;p 2 ;p 3 ;p 4 ) + T^ baba (p 1 ;p 2 ;p 3 ;p 4 ) 

+ T^ baab (p l]P2 ;p 3 ;p A ), (4.13) 

for any b ^ a. The structure of T^ 4 ) is encoded in the Dyson-Schwinger equation 
(4.14) (see also [26, 57]). In the latter the sum Y$=i schematically represents a 
summation over s, t and u scattering channels. For clarity's sake the internal indices 
are suppresses. Similarly as before, we can argue that both the third and fourth 
graphs contribute at most N~ 2 , whilst the second ("fish") graph may contribute up 
to order iV -1 . So in the large- N limit the last three diagrams may be neglected, 
provided we keep in the 4-point vertex function only graphs proportional to iV -1 . 
However, the former can be only fulfilled if we retain such a "fish" graph where 
summation over internal index on the loop is allowed. Remaining two graphs in the 

12 Actually, Ward's identities read [57] / d D x <f> b (x) = Jd D x <j> a (x) (here <f> a = fJ£; 

W is the generating functional of connected Green's functions). Performing successive variations 
with respect to (f> a {v), 4> a {z), 4> a {y) and 4>b{w), taking the Fourier transform, and setting the physical 
condition (f> c — 0, we get directly (4.13). 
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sum Yh=i (i- e -5 t and u scattering channels) are suppressed by the factor N 1 as the 
internal index on the loop is fixed. In this way we are left with the relation 



p(4) aabb^ = q) 

c+b 



d ° q bbcc (s) 
(2tt)^ [S) 



A X X r (N-l) /-i 



N 



2N 2 



L dX I (2tt 



(q 2 — m 2 + ie) ((q — Q) 2 — m 2 + ie) 
d D q i 



s=0 



(2tt) d (q 2 — m 2 + x(l — x)s + ie)' 



s=0 



(4.15) 



with Q = pi + p 2 and s — Q 2 , p±,p 2 are the external momenta. To leading order in 
1/N we may equivalently write 



A r = A + A A r M'(m 2 ) 



(4.16) 



the prime means differentiation with respect to m 2 ; Ai(m 2 ) is defined by (4.3). 
Evaluating explicitly A4'(m 2 ), we get from (4.16) 



A = 



A, 



l-A r r(2-f) (m r ) D " 4 /2 (4tt)t 



(4.17) 



Assuming that both A > and A r > (note A r < would be incompatible with 
A > and m 2 > 0), we can infer from (4.17) that 



0<A r <«^ 

- r " r(2-f) 
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(4.18) 



and so for D = 4 we inevitably get that A r = 0. The latter indicates that the theory 
is trivial [1, 25, 33], or, in other words, the O(N) <ft 4 theory is a renormalized free 
theory in the large-iV limit. This conclusion is also consistent with the observation 
that the theory does not posses any non-trivial UV fixed point in the large-iV limit 
[25, 33, 34, 60]. 

On the other hand, if we were assuming that Ao < 0, we would get indeed a 
non-trivial renormalized field theory in D = 4 (actually, from (4.17) we see that 
A — > 0_ , provided that A r is fixed and positive and D — > 4_). However, as it 
was pointed out in refs. [1, 25, 33, 35], such a theory is intrinsically unstable as the 
ground-state energy is unbounded from below. This is reflected, for instance, in the 
existence of tachyons in the theory [1, 33, 35, 61], therefore the case with negative 
Ao is clearly inconsistent. 

The straightforward remedy for this situation was suggested by Bardeen and 
Moshe [33]. They showed that the only meaningful (stable) O(N) 4 theory in the 
large- A?" limit is that with A r , Ao > 0. This is provided that we view it as an effective 
field theory at momenta scale small compared to a fixed UV cut-off A. The cut-off 
itself is further determined by (4.16) because in that case (assuming m r <C A) 

= ^ — , (4.19) 

which implies that for A r , A > we have A 2 < m 2 exp(^~). The case A 2 = 
m 2 exp( 3 y^-) corresponds to the Landau ghost [62] (tachyon pole [1, 33]). For 
reasonably small A r , A is truly huge 13 and so it does not represent any significant 
restriction. 

In passing we may observe that from (4.1) and (4.19) 

and so the fraction m 2 /A r is renormalization invariant. It was argued in [33] that 
for the part of the parameter space where A > and m 2 /A r > the ground state is 
O(N) symmetric, Goldstone phenomena cannot materialize and hence the expecta- 
tion value of the field is zero. The latter fact has been implicitly used, for instance, 
in the derivation our Dyson-Schwinger equations. To avoid a delicate discussion of 
the phase structure of the O(N) 4 theory and to emphasize our primary objective 
- i.e. hydrostatic pressure calculation, we confine ourselves to the parameter space 
defined above. Such an effective theory will provide a suitable playground to explore 

13 For example, if A r = 1 and m r w lOOMeV, we get A < 10 141 MeV or equivalently A < 10 131 K 
(this is far beyond the Planck temperature - 10 32 K). 
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all the basic salient points involved in the hydrostatic pressure calculation. Further- 
more, because the mass-shift equation (gap equation) has a particularly simple form 
in this case the hight-T analysis of the hydrostatic pressure will be easy to perform. 



4.3 Resumed hydrostatic pressure 

The partition function Z has a well known path-integral representation at finite 
temperature, namely 

Z[T] = exp($[T]) =Jv<f> exp(^[0; T]) , 

S[<j);T] = f d D xC{x). (4.21) 

Here $ = — f3Q is a Massieu function (the Legendre transform of the entropy) 
[5, 41, 42, 43] and J c d D x = J c dx J v cP _1 x with the subscript C suggesting that 
the time runs along some contour in the complex plane. In the real-time formalism, 
which we adopt throughout, the most natural version is the so called Keldysh- 
Schwinger one [5, 6], which is represented by the contour in Fig. 3. Let us mention 
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Figure 3: The Keldysh-Schwinger time path. 



that the fields within the path-integral (4.21) are further restricted by the periodic 
boundary condition (KMS condition) [5, 6, 9] which in our case reads: 

<p a {ti - i/3,x) = a (^,x) . 

As explained in Section 3, we can use for a pressure calculation the canonical energy- 
momentum tensor 0^". Employing for Q^ v (x) its explicit form (3.2) together with 
(3.6), one may write 

<0f> = f /^(2gV-^(g 2 -^))B n (g;T) + ^<T ( (E <Pl^))) , 

(4.22) 
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where E>n is the Dyson-resumed thermal propagator [5, 6], i.e. 



Note that we have exploited in (4.22) the fact that the expectation value of Q%"(x) 
is x independent. On the other hand, in (4.23) we have used the fact that m 2 . is 
q independent. In order to calculate the expectation value of the quartic term in 
Eq.(4.22), let us observe (c.f. (4.21)) that the derivative of $ with respect to the 
bare coupling Ao (taken at fixed m ) gives 



<9$[T] 

which implies that 




(4.24) 



The key point now is that we can calculate $[T] in a non-perturbative form. (The 
latter is based on the fact that we know the Dyson-resumed propagator E)n(g;T) 
(see (4.23).) Indeed, taking derivative of $ with respect to wig (keeping Ao fixed) we 
obtain 

d$[T\ iN f D I 2 v PVN f d D q 

= -(3VN M T {m 2 {T)), (4.26) 

roc 

$[T; A ; ml] = (3VN / dm 2 M T (m 2 r (T)) + $[T; A ; oo] . (4.27) 

J mi 



thus 



Let us note that 3>[T; Ao; oo] is actually zero 14 because $[T; Ao; Wq] has the standard 
loop expansion [5, 57] depicted in Fig.4. It is worth mentioning that in the latter 
expansion one must always have at least one type-1 vertex [5]. The RHS of Fig.4 
clearly tends to zero for m — > oo as all the (free) thermal propagators from which 
the individual diagrams are constructed tend to zero in this limit. The former result 
can be also deduced from the CJT effective action formalism [27] or from a heuristic 
argumentation based on a thermodynamic pressure [1]. Note that in the large-iV 
limit the fourth and fifth diagrams in Fig.4 must be omitted. 



The expectation value (4.25) can be now explicitly written as 

= 8^ 2 / dml ^^rlm L ^ —A • (4.28) 
Jmi J 2tt) d e*>P - 1 \ q 2 - ml + ie 2 v ; 




14 To be precise, we should also include in Fig.4 an (infinite) circle diagram corresponding to the 
free pressure [5, 40]. However, the later is Ao independent (although mo dependent) and so it is 
irrelevant for the successive discussion (c.f. Eq.(4.25)). 
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Figure 4: First few bubble diagrams in the $ expansion. 



In fact, the differentiation of the proper self-energy in (4.28) can be carried out 
easily. Using (4.10), we get 



+ X M' T 



8X A dX dX A (l - X M' T ) 

From Eq.(4.10) it directly follows that 

alm 2 r {T) _ 1 
alml ~ (1 - A M' T ) ' 

which, together with the definition of Mt, gives 



= - 8iV 2 f dm? M T (*P 8MtW) 



m?(T) ' * " (9m 2 . 

= 4iV 2 M 2 T {m 2 r {T)) , (4.29) 

where we have exploited in the last line the fact that .My(m 2 — > oo) = 0. Let us 
mention that the crucial point in the previous manipulations was that m r is both 
real and momentum independent. Collecting our results together, we can write for 
the hydrostatic pressure per particle (cf. Eq(3.30)) 

V(T)-V(0) = - {D ^ 1)N (fe)-(0fe|0)) 

1 (S 2 K(T))-S 2 K 2 (0))) . (4.30) 



2A 
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Applying Green's theorem to the last two integrals and eliminating the surface terms 
(for details see Appendix A) we find 

^(S|K 2 (T))-S 2 K(0))) 



2A 

= Af T (m 2 (T)) - A/"K 2 (0)) + ^ (S^K(T)) - S 2 K 2 (0))) , (4.31) 
where we have introduced new functions J\fT(m 2 (T) and A/"(m 2 ); 

W:K(T)) = I/^-£W_ %2 - m;( r)) 

A^(mJ) = Jim AfrK(T)) . (4.32) 

Eq.(4.31) can be rephrased into a form which exhibits an explicit independence of 
bar quantities. Using the trivial identity: 

2^ (^K(T))-S 2 K(0))) 

= 2^ N m *™ ~ S K 2 (°))) (z T (m 2 r (T)) + S(m 2 (0))) 

= 5 -^l (M T (ml(T)) + M(m 2 r (0))) , (4.33) 

we get 

V(T)-V(0) 

= M T {ml{T))-U{m 2 M) + 5 -^^ (M T (m 2 r (T)) + M(m 2 r (0))) , (4.34) 

where Sm 2 (T) = m 2 (T) — m 2 (0). Let us finally mention that the finding (4.30) is 
an original result of this paper. The result (4.34) has been previously obtained by 
authors [1] in the purely thermodynamic pressure framework. 



5 High— temperature expansion of the hydrostatic 
pressure in D = 4 

In order to obtain the high-temperature expansion of the pressure in D = 4, it is 
presumably the easiest to go back to equation (4.30) and employ identity (4.33). Let 
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us split this task into two parts. We firstly evaluate the integrals with potentially 
UV divergent parts using the dimensional regularization. The remaining integrals, 
with the Bose-Einstein distribution insertion, are safe of UV singularities and can 
be computed by means of the Mellin transform technique. 

Inspecting (4.30) and (4.33), we observe that the only UV divergent contributions 
come from the integrals: 

+ (tT^T) / (2^ ^ " m ' (T)) " ^ " m2(0)) ) 

5m 2 (T) r d D q ( i i \ 

4 J (2n) D \q 2 - m 2 (T) + ie q 2 -m 2 (0)+^J ' { ' } 
which, if integrated over, give 

(5.i)= + r(^)r(f + 1) , o _ o v 

( J D-l)r(^i)(47r)T VV rV > 



Sm 2 (T) T(l - f ) 



+ 2 



2 



((^( T )) f_1 + (^r(0)) f_1 ) • (5.2) 



4(4tt) 

Taking the limit D = 4 — 2e — > 4 and using expansions 

a x+£ = a x (l + e lna + 0(e 2 )) , 
(7 is the Euler-Mascheroni constant) we are finally left with 

The fact that we get finite result should not be surprising as entire analysis of Section 
3 was made to show that V(T) — V(0) defined via Q^ v is finite in D = 4. 

We may now concentrate on the remaining terms in (4.30), the latter read (we 
might, and we shall, from now on work in D = 4) 

(5.4) 

Our following strategy is based on the observation that the previous integrals have 
generic form: 



rJ 1 



(2tt) 3 1 el^ol/3 _ 1 

— / dx(x 2 -l)— -, 5.5 

7T 2 Ji K ' e*y - 1 v ' 
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with v — 0, 1 and y = m r j3. Unfortunately, the integral (5.5) cannot be evaluated 
exactly, however, its small y (i.e. high-temperature) behavior can be successfully 
analyzed by means of the Mellin transform technique [5, 40]. Before going further, 
let us briefly outline the basic steps needed for such a small y expansion. 

The Mellin transform f(s) is done by the prescription [5, 40, 63, 64, 65, 66]: 

roo 

f(s) = / dx X s - 1 f(x) , (5.6) 
Jo 

with s being a complex number. One can easily check that the inverse Mellin 
transform reads 

I rioo+a 

/(*) = -7^-J ds X ~ S > ( 5 - 7 ) 

tyZlTj J-ioo+a 

where the real constant a is chosen in such a way that f(s) is convergent in the 



l 



e x v-l 



neighborhood of a straight line (— too + a, too + a). So particularly if f(x) 
one can find ( [65]; formula 1.3.19) that 

f(s) = r(s){(s)y-, (Res>l), (5.8) 

where ( is the Riemann zeta function (((s) = J2^=i n~ s ). Now we insert the Mellin 
transform of f(x) = eX ]_ l to (5.5) and interchange integrals (this is legitimate only 
if the integrals are convergent before the interchange). As a result we have 

roo I pioo+a /jo 

/ dx g(x) -— = / -— T(sK(s)y' s g(l - s) , (5.9) 

JO e xy — 1 J-ioo+a l{ln) 

with g(x) = 9(x — 1) [x 2 — l) 1 ^. Using the tabulated result ( [66]; formula 6.2.32) 
we find 

g(l - s) = \b{-v - 1 + \s- § + v) , (Res>2 + 2z/), (5.10) 

with B( ; ) being the beta function. Because the integrand on the RHS of (5.9) is 
analytic for Res > 2 + 2v and the LHS is finite, we must choose such a that the 
integration is defined. The foregoing is achieved choosing a > 2 + 2v. Another useful 
expressions for g(l — s) are ( [66]; formula 1.2.34 or 1.2.37) 

g{l-s) = B(l + u;-2-2u + s) 2 F 1 [-\ - v; -2 - 2v + s; -\ - v + s; -1] 
= 2^ +u B(\ + v--2-2u + s) 2 F 1 [-\ -v;\ + v;-\-v + s;\], 

where 2 F\ is the (Gauss) hypergeometric function [66]. Using identity 

2 2x ~ 1 

r(2x) = -j=- r(x)r{x + \) , 
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we can write 



(5.9) = i^±^ r * -^n\s)as)(\ y y S T{-v-l + \s). (5.11) 

4^/7T J-joo+a l{2lT) v 7 

The integrand of (5.11) has simple poles in s = —2n {n — 1, 2, . . .), s — 1, s — 
—2n + 2v + 2 (n — 0, 1, . . . , v) and double pole in s = 0. An important point in the 
former pole analysis was the fact that £(s) has simple zeros in —2m (m > 0) and 
only one simple pole in s — 1. The former together with identity 

_ ( X \ x u . . „ / 1 — X\ x — l 

r 77 ^(W=r — 7T— c(i-x), 



shows that no double pole except for s = is present in (5.11). Now, we can close 
the contour to the left as the value of the contour integral around the large arc is zero 
in the limit of infinite radius (c.f. [65] and [67]; formula 8.328.1). Using successively 
the Cauchy theorem we obtain 

V? (5.9) 

^ an-2,-2 *- 2n+2 » +2 {-n + v)\{-lY\B_ 2n+2v+2 \ 
^ Q y n\ {-In + 2v + 2)! 2 4 "- 4 ^ 4 

~ 2TO vr- 2 " (2n)!C(l + 2n) (-1)"+^ 
n! (n + l + z/)!2 4 «-! 

I ,-i ^ Llg + ^ 2 " +3 | 2 Llg [in r M | -,- i V 1 1 (512) 
+ y (2^ + 2)! + (,/ + !)! V n UJ +7 2 hk} ,{ } 

where B^s are the Bernoulli numbers. Let us mention that for ((2n + 1) only 
numerical values are available. 

Inserting (5.12) back to (5.4), we get for V(T) - V(0) 

V(T) - P(0) = (5.3) + i h{m r {T)) + h{m r {T)) 

_T^_T^( 2(rr , 6m?(T) \ Tm r (T) (m*(T) 6m\T) \ 
~ 90 24 \ mAI) 2 ) + 4tt v 3 4 ) 

m 2 (T)m 2 (0) ( (m r (0)\ l\ 
+ 32 k* { ln {-T^) +1 -2)-128^ 

- V fm 2 m - fe±2 5m 2 m^ m^(T) ^ (2n)! C(l + 2n) (-l)H- 
^ A J 2 1 " T 2 «n! (n + 2)!2 4 ™+ 4 

(5.13) 
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Note that (5.3) cancelled against the same term in | I 2 (m r (T)) + — ^ I (m r (T)). 
One can see that (5.13) rapidly converges for large T, so that only first four terms 
dominate at sufficiently high temperature. The aforementioned terms come from the 
poles nearby the straight line (— ioo+a, ioo+a) (the more dominant contribution the 
closer pole). It is a typical feature of the Mellin transform technique that integrals 
of type 

dx g{x) 



Jo 




jo e*y - 1 ' 

can be expressed as an expansion which rapidly converges for small y (high-temperature 
expansion) or large y (low-temperature expansion) 15 . 

Expansion (5.13) is the sought result. To check its consistency we will apply it 
to two important cases: high T case and m r (0) = case. Concerning the first case, 
note that for a sufficiently large T we can use the high-temperature expansion of 
8m 2 (T) found in Appendix B. Inserting (B.6) to (5.13) we obtain 

V(T) - V(0) = — - T2mliT) + T * mriT) 
y ' y ' 90 24 12vr 

gMQ + m2 f>) + Q U T) In f^|) . (5.14) 
24tt 16tt 2 J y \ TAn ) ) y ' 

Up to a sign, the result (5.14) coincides with that found by Amelino-Camelia and 

Pi [2] for the effective potential 16 . Actually, they used instead of the N — > oo limit 
the Hartree-Fock approximation which is supposed to give the same V e ff as the 
leading 1/N approximation [62]. 

As for the second case, we may observe that our discussion of the mass renor- 
malization in Section 3.1 can be directly extended to the case when m r (0) = (this 
does not apply to our discussion of A r !). Latter can be also seen from the fact that 
(5.13) is continuous in m r (0) = (however not analytic). The foregoing implies that 
the original massless scalar particles acquire the thermal mass m 2 (T) = 5m 2 (T) . 
From (5.13) one then may immediately deduce the pressure for massless fields a in 
terms of Sm(T). The latter reads 

V (T)-V(0) = ^ - ^ {6m{T))2 I T (MT))3 



90 48 48 7i 

~ (5m(T)) 2 ^ 7T- 2 - 2 (2n)! C(l + gn) (-1)" +1 
2^ T 2n (n — 1)! (n + 2)! 2 4n+5 ' 1 J 



15 By the same token we get the low-temperature expansion if the integral contour must be closed 
to the right. 

16 Let us remind [2, 27, 31] that from the definition of V e ff the thermodynamic pressure is —V e ff. 
In order to obtain (5.14) from V e ff in [2], one must subtract the zero temperature value of V e /f 
and restrict oneself to vanishing field expectation value and positive bare mass squared. 
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This result is identical to that found by Drummond et al. in [1]. 

A noteworthy observation is that when the energy of a thermal motion is much 
higher then the mass of particles in the rest, then the massive theory approaches 
the massless one. This is justified in the first (high-temperature dominant) term 
of (5.13) and (5.15). This term is nothing but a half of the black body radiation 
pressure for photons [41, 42] (photons have two degrees of freedom connected with 
two transverse polarizations). One could also obtain the temperature dominant 
contributions directly from the Stefan-Boltzmann law [5, 41, 42] for the density 
energy (i.e. (O 00 )). The formal argument leading to this statement is based on the 
noticing that at high energy (temperature) the scalar field theory is approximately 
conformally invariant, which in turn implies that the energy-momentum tensor is 
traceless [52]. Taking into account the definition of the hydrostatic pressure (2.8), 
we can with a little effort recover the leading high-temperature contributions for 
the massive case. 

6 Conclusions 

In the present article we have clarified the status of hydrostatic pressure in (equi- 
librium) thermal QFT. The former is explained in terms of the thermal expectation 
value of the "weighted" space-like trace of the energy-momentum tensor 6^. In 
classical field theory there is a clear microscopic picture of the hydrostatic pressure 
which is further enhanced by a mathematical connection (through the virial theo- 
rem) with the thermodynamic pressure. In addition, it is the hydrostatic pressure 
which can be naturally extended to a non-equilibrium medium. Quantum theoretic 
treatment of the hydrostatic pressure is however pretty delicate. In order to get a 
sensible, finite answer we must give up the idea of total hydrostatic pressure. In- 
stead, thermal interaction pressure or/and interaction pressure must be used (see 
(3.30) and (3.31)). We have established this result for a special case when the the- 
ory in question is the scalar 4 theory with O(N) internal symmetry; but it can 
be easily extended to more complex situations. Moreover, due to a lucky interplay 
between the conservation of 6 Miy and the space-time translational invariance of an 
equilibrium (and T = 0) expectation value we can use the simple canonical (i.e. 
unrenormalized) energy-momentum tensor. In the course of our treatment in Sec- 
tion 3 we heavily relied on the counterterm renormalization, which seems to be the 
most natural when one discusses renormalization of composite Green's functions. To 
be specific, we have resorted to the minimal subtraction scheme which has proved 
useful in several technical points. 

We have applied the prescriptions obtained for the QFT hydrostatic pressure to 
4 theory in the-large N limit. The former has the undeniable advantage of being 
exactly soluble. This is because of the fact that the large- N limit eliminates "nasty" 
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classes of diagrams in the thermal Dyson-Schwinger expansion. The survived class 
of diagrams (superdaisy diagrams) can be exactly resumed, because the (thermal) 
proper self-energy E, as well as the renormalized coupling constant A r are momen- 
tum independent. We have also stressed that the O(N) 4 theory in the large-iV 
limit is consistent only if we view it as an effective field theory. Fortunately, the up- 
per bound on the UV cut-off is truly huge, and it does not represent any significant 
restriction. For the model at hand the resumed form of the pressure with m r (0) = 
was firstly derived (in the purely thermodynamic pressure context) by Drummond 
et al. in [1]. We have checked, using the prescription (3.30) for the thermal inter- 
action pressure, that their results are in agreement with ours. The former is a nice 
vindication of the validity of the virial theorem for the QFT system at hand. In 
this connection we should perhaps mention that the latter is by no means obvious. 
For example, for quantized gauge fields the conformal (trace) anomaly may even 
invalidate the virial theorem [5]. The fact that this point is indeed non-trivial is 
illustrated on the QCD case in [68]. 

The expression for the pressure obtained was in a suitable form which allowed 
us to take advantage of the Mellin transform technique. We were then able to write 
down the high-temperature expansion for the pressure in D = 4 (both for massive 
and massless fields) in terms of renormalized masses m r (T) and m r (0). We have 
explicitly checked that all UV divergences present in the individual thermal diagrams 
"miraculously" cancel in accordance with our analysis of the composite operators in 
Section 3. 
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A Appendix 

In this Appendix we give some details of the derivation of Eq.(4.31). We particularly 
show that the surface integrals arisen during the transition from (4.30) to (4.31) 
mutually cancel among themselves. As usual, the integrals will be evaluated for 
integer values of D and corresponding results then analytically continued to a desired 
(generally complex) D. 

The key quantity in question is 



+ 
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- 5/(2^1 ((Z^I))^ 2 -^ ')' < A '« 

Applying Green's theorem (i.e. integrating by parts with respect to q) on (A.l) one 
finds 

(A.l) = JVtW(T)) -M(m 2 J0)) 

- & ^by / L s -= ds 11 "<« 2 - m2(0 » »<*> ■ (A - 2) 



As usual, ab = Y,f=i a «bj and S R 2 is a (Z? — 2)-sphere with the radius i?. The 
expressions for Af T and A/" are done by (4.32). 

With the relation (A. 3) we can show that the surface terms cancel in the large 
R limit. Let us first observe that 

d l0 f j n, 2 2/m\\ 2 %o) 



lim / - q ° n 1 / dsq% 2 -m?(T)) . 

27rVig D - 1 f_dqo_ af2 p2 2m , 2% ) 



7T 2 ^ .... 

hm -. r- / , dan —3 = 0. (A. 3) 



R^oo 



In 2-nd line we have exploited Gauss's theorem and in the last line we have used 
L'Hopital's rule as the expression is in the indeterminate form 0/0. The remaining 
surface terms in (A. 3) read 

ii m J ds q { 9{q2 - m ' (T)) - 9{q2 - m ' (o)) ) 

= lim =f—^- J / , - / , \dq = 0. (A.4) 

R^OO 2 D-2 T (D_1\ \J^ R 2 +m 2 [T) J^ R 2 +m 2 {0) J 

The last identity follows either by applying L'Hospital's rule or by a simple transfor- 
mation of variables which renders both integrals inside of {. . .} equal. Expressions 
on the last lines in (A. 3) and (A.4) can be clearly (single-valuedly) continued to the 
region ReD > 1 as they are analytic there. We thus end up with the statement that 

(A.l)=Ar T K 2 (T))-ATK 2 (0)). 
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B Appendix 

In this appendix we shall derive the high-temperature expansion of the mass shift 
Sm 2 (T) in the case when fields a are massive (i.e. mf(0) ^ 0). 

Consider Eqs.(4.1) and (4.10). If we combine them together, we get easily the 
following transcendental equation for Sm 2 (T) 

5m 2 (T) = A |^K(T)) - M(m 2 (0)) + h (m 2 (0) + <5m 2 (T))j . (B.l) 

Here M and I are done by (4.3) and (5.5), respectively. 

Now, both Ao and M. are divergent in D = 4. If we reexpress Ao in terms of A r , 
divergences must cancel, as 5m 2 (T) is finite in D = 4. The latter can be easily seen 
if we Taylor expand M, i.e. 

M(m 2 (T)) = M(m 2 (0)) + 5m 2 {T) M'(m 2 (0)) + M(m 2 r (0); 5m 2 {T)) . (B.2) 

Obviously, Xi is finite in D = 4 as M. is quadratically divergent. Inserting (B.2) to 
(B.l) and employing Eq.(4.16) we get 

5m 2 (T) = A r {-M(m 2 (0); 5m 2 (T)) + ^/ K 2 (0) + 5m 2 (T))| . (B.3) 

This is sometimes referred to as the renormalized gap equation. In order to deter- 
mine Ai we must go back to (B.2). From the former we read off that 

M(m 2 (T);5m 2 (T)) 



M(m 2 (T)) - M(m 2 r (0)) - 5m 2 (T) M'(m 2 (0)) 

' (1 ~$- (K(T))^ 1 - (m 2 (0))f - 1 - 6m 2 (T)£ - 1) (m 2 (0))f-2 



32,n m ' (r)ln U(5) " W(r) ' (R4) 



So 



m 2 (0)+<5m 2 (T)) lnfl + ^S) -5m 2 (T) i 
^ 2 (T) = A r < — 12 m ^ I 21 + (B.5) 

Analogous relation was also derived in [32] where authors used finite temperature 
renormalization group. In the latter the zero-momentum renormalization prescrip- 
tion was utilized. Eq.(B.5) was firstly obtained and numerically solved in [1]. It was 
shown that the solution is double valued. The former behavior was also observed 
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in the effective action approach. Namely by Abbott et al. [35] at T = 0, and by 
Bardeen and Moshe [33] at both T = and T ^ 0. The relevant solution is only 
that which fulfils the consistency condition 5m 2 (T) — > when T — > 0. For such a 

solution it can be shown (c.f. [1], Fig. 3) that S ™2(Q) <S 1 for a sufficiently high T. 
So the high-temperature expansion of (B.5) reads 



Sm 2 (T) = \ r 



(<5m 2 (T)) 2 _ (<5rn 2 (T)) 3 (<5m 2 (T)) 
2m 2 (0) Qmf(O) 12mp(0) 

32^ 



4 




= = ^ - A ^ (T) T + Q U(T) In f ^1 . (B.6) 
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